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Abstract 

In this article, we propose the application of a modified Taylor series method (MTSM) for the approximation of 
nonlinear problems described on finite intervals. The issue of Taylor series method with mixed boundary conditions is 
circumvented using shooting constants and extra derivatives of the problem. In order to show the benefits of this 
proposal, three different kinds of problems are solved: three-point boundary valued problem (BVP) of third-order with 
a hyperbolic sine nonlinearity, two-point BVP for a second-order nonlinear differential equation with an exponential 
nonlinearity, and a two-point BVP for a third-order nonlinear differential equation with a radical nonlinearity. The result 
shows that the MTSM method is capable to generate easily computable and highly accurate approximations for 
nonlinear equations. 
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Introduction 

Nonlinear differential equations are a popular tool to 
model complicated dynamical phenomenons of many 
branches of sciences. Unfortunately, a drawback of this 
approach arises when exact solutions are required but 
not available. On one side, pure numerical methods 
employed to solve nonlinear differential equations can 
exhibit numerical instabilities, oscillations or false equi- 
librium states, among others. On the other side, approx- 
imative methods are a good option when semi-analytic 
solutions are required. Some examples of such methods 
are: homotopy perturbation method (HPM) (Filobello- 
Nino et al. 2012a; Filobello-Nino et al. 2012b; He 1999; 
2009; Khan et al. 2013; Vazquez-Leal 2012; Vazquez-Leal 
et al. 2012a; Vazquez-Leal et al. 2012b; Vazquez-Leal et 
al. 2012c), homotopy analysis method (HAM) (Hassana 
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and El-Tawil 2011; He 2004; Tan and Abbasbandy 2008), 
variational iteration method (VIM) (Chang 2010; Khan et 
al. 2012), Taylor series method (TSM) (Barrio et al. 2011; 
Rodriguez and Barrio 2012; Shiraishi et al. 2011; Wazwaz 
1998), Adomian decomposition method (Duan and Rach 
2011; Wazwaz 1998), among others. Nevertheless, TSM 
highlights because of its simplicity and power; it does 
not require a perturbation parameter as the perturba- 
tion based techniques or trial functions as HAM or HPM 
does. In addition, TSM is straightforward and can be pro- 
grammed using computer algebra packages like Maple 
or Mathematica. What is more, TSM method was con- 
ceived as a tool to solve differential equations governed 
by Dirichlet conditions (DC), although, mixed boundary 
condition (MBC) problems are important and common in 
several fields of physics. Therefore, in this work, we pro- 
pose a modified Taylor series method (MTSM) for MBC 
problems that is based on the induction of shooting con- 
stants (SC) (Stoer and Bulirsch 2002). Such constants arise 
from two sources: 



Springer 



© 2014 Vazquez-Leal et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons 
Attribution License (http://creativecommons.0rg/licenses/by/2.O}, which permits unrestricted use, distribution, and reproduction 
in any medium, provided the original work is properly credited. 



Vazquez-Leal etal. SpringerPlus 2014, 3:1 60 
http://www.springerplus.eom/content/3/1 /1 60 



Page 2 of 7 



1. The conversion of MBC to DC. This process implies 
the choosing of an expansion point (usually at zero) 
and the conversion of the MBC that are not at the 
expansion point to SC constants to be determined 
later by MTSM. The number of these constants is 
the same as the number of boundary conditions that 
are not at the expansion point. 

2. Increasing the order of the original differential 
equation by the application of extra derivatives; as a 
strategy to add integration constants to solution that 
work as shooting constants or adjustment 
parameters. The number of this SC constants 
depends of the level of accuracy that we require from 
the approximate MTSM solution. For all the cases 
study of this work only one extra derivative is 
required to obtain a good fitting with respect to the 
exact solution, although it depends of the particular 
problem under study. 

The aforementioned combined shooting technique of 
MTSM aids to circumvent the issue of TSM method 
with MBC. In order to show the benefits of this pro- 
posal, three nonlinear problems described with MBC on 
finite intervals are solved: three-point BVP for a third- 
order nonlinear differential equation with a hyperbolic 
sine nonlinearity (Duan and Rach 2011), two-point BVP 
for a second-order nonlinear differential equation with 
an exponential nonlinearity (Duan and Rach 2011; Scott 
and Vandevender 1975) and a two-point BVP for the 
third-order nonlinear differential equation with a radical 
nonlinearity (Duan and Rach 2011). 

This paper is organized as follows. In Section 'MTSM 
method', we introduce the basic idea of MTSM method. 
In Section 'Cases study', we show the solution proce- 
dure for three nonlinear problems. Numerical simula- 
tions and a discussion about the results are provided in 
Section 'Numerical simulation and discussion'. Finally, a 
brief conclusion is given in Section 'Conclusion! 

MTSM method 

We consider a nonlinear differential equation expressed as 



Firstly, we increase the order of the differential equation 



d k 

u (n+k) = [Ar(M) _f( x )] t x eQ., (3) 

dx k 



where k is a constant related to the number of the desired 
SC constants. 

Next, it is possible to express the Taylor series solution 
for (3) as 



u'(X 0 ) j u"(X 0 ) 2 
Ut = U{Xq)-\ — (X — Xo) H (x — Xo) 



V. 



2! 



H — (x-xq) H — (x-x 0 ) +• 



3! 



4! 



(4) 



where Xq is the expansion point and derivatives u^{xq) 
(i = 0, 1, . . . ) are expressed in terms of the parameters and 
boundary conditions of (3). 

As we require to solve MBC problems, the boundary 
conditions not located at the chosen expansion point Xo 
will be replaced by shooting constants giving as result tra- 
ditional DC conditions. Next, in order to obtain the coef- 
ficients of (4) (u (i) (xo), i = 0,1,...), MTSM requires (I) 
calculate the successive derivatives of (3) and (II) evaluate 
each derivative using the Dirichlet conditions. Finally, in 
order to fulfil the boundary conditions originally replaced 
by the SC constants is necessary to evaluate (4) in such 
points; then, the resulting system of equations is solved 
to obtain the value of the SC constants. It is important 
to remark that the order of the Taylor expansion (4) is 
chosen in order to include all the shooting constants in 
the polynomial; as long as we satisfy such condition the 
order of the Taylor expansion can be increased to improve 
accuracy. 

The constants due to the extra /c-derivatives (see (3)) are 
applied to minimize the mean square residual (MSR) error 
defined as 



u (n) =N{u)-f{x), xeQ., (1) J \ u t -N(u T )+f(x)j dx, 



(5) 



having as boundary condition 

x e r, 



du 

B(u, — ) = 0, 
dr] 



(2) 



where n is the order of the differential equation, AT is a 
general operator; f(x) is a known analytic function, B is a 
boundary operator, Y is the boundary of domain Q., and 
du/dr] denotes differentiation along the normal drawn 
outwards from Q. 



where ut is the approximated TSM solution (4), and 
[Xi,Xf] is the finite interval delimited by the MBC. 

Cases study 

In the present section, we will solve three cases study to 
show the utility of the MTSM method to solve nonlinear 
problems. For all cases study the expansion point of TSM 
is at #o = 0 and the derivatives were performed using 
Maple 17 software. 
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Third-order nonlinear equation 

Consider the three-point BVP for the third-order nonlin- 
ear differential equation with a hyperbolic sine nonlinear- 
ity (Duan and Rach 2011) 

u!" - 1 -*sinh(«) = 0, u(0) = 0,w(0.25) = 1,«(1) = 0, 
0 < x < 1, 

(6) 

where prime denotes derivative with respect to x and the 
exact solution is unknown. 
Firstly, we derive (6), resulting 

u>' v ' — xu' cosh(w) — sinh(w) = 0, 

u(0) = 0,u'(0) = c 1: u"(0) = c 2 ,u"'(0) = c 3 ,0 < x < 1, 

(7) 

where the boundary conditions of (6) are replaced by its 
Dirichlet conditions accordingly to the increased order. 

Next, we resolve (7) for m ( ' v) and perform successive 
derivatives, resulting 

w ( ' v) = xu' cosh(w) + sinh(w), 

w (v) = 2w' cosh(w) + x(u') 2 sinh(M) + xu" cosh(w), 

u (vi) = 3(w') 2 sinh(w) + 3u" coshO) + • • • 

(8) 



Next, c 3 is used as adjustment parameter to minimize 
the mean square residual (MSR) error by resolving 



3 

9c 3 



/ (m™ — 1 — x sinhOr)) 2 dx 
Jo 



(12) 



where c\ and C2 were previously substituted by (11). It 
is important to notice that in order to obtain a symbolic 
expression for C3 the hyperbolic sine was replaced by its 
fifth-order Taylor series. 

The result of solving (12) is c 3 = 1.1353380202 giving a 
minimum MSR error of 



/ (u% -l-x sinhOr)) 2 dx = 0.007337036421. 
Jo 



(13) 



Second-order nonlinear differential equation 

Consider the two-point BVP second-order nonlinear dif- 
ferential equation with an exponential nonlinearity (Duan 
and Rach 2011; Scott and Vandevender 1975) 

u" - exp(w) = 0, w(0) = 0, w(l) = 0, 0 < x < 1, 

(14) 

where prime denotes derivative with respect to x and 
exact solution is 



Now, the boundary conditions of (7) are substituted into 
(8), yielding 



u iiv \0) = 0, 
u (v) (0) = 2ci, 
w (W) (0) = 3c 2 , 



(9) 



Finally, using (9) and the initial conditions of (7), we can 
formulate the sixth-order Taylor series expansion (see (4)) 



1 9 1 , 1 , 1 , 
ut(x) = c\x + -C2X H — C3X H c\x H C2X . 

2 6 60 240 



(10) 



Finally, if we substitute the boundary conditions 
w(0.25) = 1 and u(l) = 0 into (10) and solve for the 
shooting constants (ci and C2), it results 



ci = 5.347659477 + 0.04132001785c 3 , 
c 2 = -10.78371002 - 0.4139015236c 3 . 



(11) 



u(x) = 21n^Csec^ C(2 * 4 1) ^-ln(2), 
C = 1.336055694906108. 



Now, we derive (14), resulting 

u'" - u' exp(w) = 0, u(0) = 0, u'(0) = c 1: u"(0) = c 2 , 
0 < x < 1, 

(16) 

where the boundary conditions of (14) are replaced by its 
Dirichlet conditions accordingly to the increased differen- 
tial equation order. 

As the aforementioned procedure for first case study, we 
obtain the following sixth-order Taylor series 



ut(x) =c\x + ]-c 2 x 2 + \c\x 3 + -^-(c\ + c 2 )* 4 
2 6 24 



1 3 s 
+ — (q + 3cic 2 + ciK 

+ + 6c i c 2 + 5cj + 3cj + c 2 )x 6 . 



(17) 
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Finally, if we substitute the boundary condition w(l) = 0 
into (17) and solve for the shooting constant C2, it results 

2 „ 391 

C2 = -ci - 3ci — 

6 

+ ~V24c7 + 144ci 3 + 4596ci 2 + 3924ci + 152881, 
6 

(18) 

where the negative square root term was discarded 
because it did not minimize the mean square residual 
error. 

Repeating the procedure for first case study (see (12)), 
we find that C\ = -0.4582864419 minimize the MSR error, 
resulting 

pi 

/ (*4 - exp( Mr )) 2 dx = 0.0004389212651, (19) 
Jo 

where ut corresponds to (17). 



Third-order nonlinear differential equation with a radical 
nonlinearity 

Consider the two-point BVP for the third-order nonlinear 
differential equation with a radical nonlinearity (Duan and 
Rach 2011) 



Firstly, we derive (20), resulting 



uf" + Vl - u 2 = 0, u(0) = 0,«'(0) = l,u{n/2) = 1 
0 < x < 7T/2, 

(20) 

where prime denotes derivative with respect to x and 
exact solution is 



0, 



m(0) = 0,«'(0) = l,u"(0) = ci,m"'(0) = c 2 , 0<x< Jt/2, 

(22) 

where the boundary conditions of (20) are replaced by its 
Dirichlet conditions accordingly to the increased order. 

As the aforementioned procedure for first case study, we 
obtain the following sixth-order Taylor series 

1,1, 1 , 1 , 

ut(x) = x H — c\x H — cqx H x H c\x . 

2 6 120 240 

(23) 

Finally, if we substitute the boundary condition 
m(7t/2) = 1 into (23) and solve for the shooting constant 
C2, it results 



C2 = 



1 ci?T 6 + 4tt 5 + 1920ci?r 2 + 76807T - 15360 



320 



(24) 



Repeating the procedure for first case study (see (12)), 
we find that c\ = 0.001065300514 minimizes the MSR 
error, resulting 



' {ri? + y 1 - u\\ dx = 0.0001811801833, 



(25) 



u(x) = sin(^). 



(21) where u t corresponds to (23). 
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Figure 2 Exact solution (1 5) (solid circles) and approximate MTSM solution (17) (solid line) for (14). The MSR error is 0.000438921 2651 . 



Numerical simulation and discussion 

From Figures 1, 2 and 3, we observe the high accu- 
racy of the proposed MTSM approximations for all cases 
study. For all cases study, only one extra derivative was 
required to obtain an acceptable low MSR error. It is 
important to mention that for Figures 2 and 3 the com- 
parison is of MTSM approximations are versus the exact 
solutions (see (15) and (21)) depicting the accuracy of the 
approximations. Nonetheless, the problem (6) do not pos- 
sess a known exact solution and the build-in professional 
numerical routines for BVP problems from Maple does 



failed to deal with this problem. Therefore, the only ref- 
erence to know the error from exact solution is the MSR 
error of 0.007337036421 (see (13)), which is considered 
very low. This means that (10) is highly accurate. 

The usefulness of coupling of a shooting method (Stoer 
and Bulirsch 2002) along with extra derivatives and the 
TSM method was exhibited by the solution of different 
highly nonlinear boundary value problems expressed in 
terms of nonlinearities such as: high order derivatives 
combined with hyperbolic sine, exponential and radi- 
cal terms, among others. What is more, the shooting 



— MTSM • Exact 




x 

Figure 3 Exact solution (21 ) (solid circles) and approximate MTSM solution (23) (solid line) of (20). The MSR error is 0.0001 81 1 801 833. 
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constants were used to fulfil the boundary conditions 
originally discarded by the artificial Dirichlet conditions. 
Finally, an extra derivative induced an extra shooting con- 
stant to minimize the MSR error giving as result high 
accurate (see (13), (19) and (25)) handy power series 
solutions. Finally, if users require more accurate approx- 
imated solutions, they should augment the number of 
SC constants (increasing k) to improve the potential of 
minimizing the MSR error (5). 

In this work, we presented a modified Taylor series 
method to deal with nonlinear problems exhibiting mixed 
boundary conditions defined on finite intervals. The 
aforementioned procedure and results show that MTSM 
can obtain power series solutions using only deriva- 
tives without requiring to solve a system of differential 
equations or the proposal of trial functions as HPM (He 
1999; 2009) or HAM (He 2004; Tan and Abbasbandy 
2008) methods, or an iterative solution procedure of inte- 
grals as VIM (Chang 2010) method. In addition, MTSM 
is not based on the existence of a perturbation parame- 
ter (Filobello-Nino et al. 2013). Therefore, further work 
will address more potential applications of the proposed 
method to other type of problems or inclusive other type 
of boundary conditions as: Robin or Neumann. 

Conclusion 

This work introduced the application of a modified Taylor 
series method (MTSM) for solving boundary value prob- 
lems (BVPs) with mixed boundary conditions defined on 
a finite interval. We were able to obtain accurate, easy 
computable, handy approximations for all cases study. The 
shooting constants arising from the substitution of the 
mixed boundary conditions by Dirichlet conditions and 
the extra derivatives of the differential equation demon- 
strate - with examples - to be a powerful strategy that 
provides easy computable and accurate approximations. 
In addition, more extra derivatives can be applied to the 
differential equation to increase the number of shoot- 
ing/adjustment constants, giving as result an enhanced 
convergence of the MTSM method. 
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